NASA TECHNICAL NOTE 


04 

On 

■ 

CD 


<C 

co 



uj 

tr 1 

tr s^S 

LOAN COPY: RETUF^^S 
AFWL (WLIL-2; S 
KIRTLAND AFB, N , 


BUCKLING OF AN EQUATORIAL 
SEGMENT OF A SPHERICAL SHELL 
LOADED BY ITS OWN WEIGHT 


by Robert E. Blum and Harvey G. McComb, Jr. 

Langley Research Center 
Langley Station , Hampton, Va. 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION • WASHINGTON, D. C. • DECEMBER 1968 


tech library kafb, nm 




TECH LIBRARY KAFB, NM 



Q131b73 
NASA TN D-4921 


BUCKLING OF AN EQUATORIAL SEGMENT OF A SPHERICAL SHELL 

LOADED BY ITS OWN WEIGHT 

By Robert E. Blum and Harvey G. McComb, Jr. 

Langley Research Center 
Langley Station, Hampton, Va. 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 


For sale by the Clearinghouse for Federal Scientific and Technical Information 
Springfield, Virginia 22151 — CFSTI price $3.00 



BUCKLING OF AN EQUATORIAL SEGMENT OF A SPHERICAL SHELL 
LOADED BY ITS OWN WEIGHT 

By Robert E. Blum and Harvey G. McComb, Jr. 

Langley Research Center 

SUMMARY 

Buckling equations are solved for a simply supported equatorial segment of a 
spherical shell loaded in the axial direction by its own weight. The Galerkin method is 
used to obtain a critical-thickness parameter for buckling as a function of the shell mate- 
rial and geometric properties and of the relative proportions of the shell weight supported 
at the upper and lower edges. Results of the present method are compared with those of 
a simplified approach in which the spherical shell segment is assumed to buckle if the 
maximum compressive stress is greater than the critical compressive stress for a com- 
plete spherical shell loaded by uniform external pressure. This comparison indicates 
that the simplified approach yields larger thicknesses than the method presented in the 
present report. The particular proportions of the shell weight which should be supported 
at the upper and lower edges to minimize the tendency to buckle are also determined. 
Results have been used in the design of a large rectangular reinforced-plastic shell seg- 
ment which simulates the lunar surface in the lunar orbit and landing approach simulator 
at the Langley Research Center. 


INTRODUCTION 

The lunar orbit and landing approach simulator at the Langley Research Center 
incorporates a large rectangular segment of a thin-walled spherical shell which simulates 
the lunar surface. (See fig. 1.) This shell segment is loaded by its own weight and can be 
supported completely at either its upper or lower edge or can be supported in various 
proportions at these edges. (See fig. 2.) To design this shell segment properly, its sta- 
bility must be considered. The stability of the rectangular segment can be approximated 
by the solution to the buckling equations for a complete circumferential segment of a 
spherical shell centered at the equator (equatorial segment) and loaded by its own weight; 
this problem is considered herein. 

A simplified approach to the problem is to assume that the spherical shell segment 
buckles if the maximum compressive stress, as predicted by linear membrane theory, is 
greater than the critical compressive stress for a complete spherical shell loaded by 



uniform external pressure. The purpose of this paper is to present a more realistic 
solution, to compare this solution with the simplified approach, and to determine how best 
to distribute the supporting line reactions at the upper and lower edges of the shell so that 
buckling tendencies are minimized. 

In the present paper a modified form of the buckling equations for shallow shells as 
given in reference 1 is utilized. These basic equations can be derived by a method given 
in reference 2 in which the shallow shell is represented as a flat plate with an initial 
deflection. To solve the present problem, the shell geometry is specialized to an equa- 
torial segment of a spherical shell, and the Galerkin method is used to solve the resulting 
equations for a critical-thickness parameter. Plots are given showing the critical- 
thickness parameter for buckling as a function of shell material and geometric properties 
and of support conditions and comparing the present results with results of the simplified 
approach mentioned previously. 


SYMBOLS 

A m ,Aq,A M constants . .,Am) 

{a} ,{b}> (M x 1) column matrices 

pqm) »[}>qm] X matrices 

C constant of integration 

Et^ 

D flexural stiffness of shell wall, 

12(1 -m 2 ) 

E Young’s modulus 

I unity matrix 

k proportion of weight of equatorial segment supported at upper edge 

L length of equatorial segment 

M number of terms in assumed truncated series solution 

m,q integers 
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N 


integer 


Nx,Ny,Nxy 

n 

Px’Py’Pn 

R 

t 

u,v,w 

x,y,z 

a 

j3 = R/L 
r 

V 2 ,V 4 

X 

* L cr 
X* 

M 

$ 


membrane stress resultants (see fig. 3) 

number of complete circumferential waves 

components of body force acting in the meridional (also in the axial), 
circumferential, and normal directions, respectively 

radius of equatorial segment 

shell thickness 

meridional, circumferential, and normal displacements of the spherical shell, 
respectively 

axial, circumferential, and radial coordinates, respectively 
specific -weight parameter, yL/E 


weight per unit volume of shell-wall material (on earth, gravitational 
acceleration is 32.17 ft/sec 2 or 9.807 m/sec 2 ) 

Laplacian and bi-Laplacian operators, respectively 


thickness parameter, 


t* 


12(1 - jli 2 )r 2 

critical-thickness parameter obtained by the present method 
critical-thickness parameter obtained by the simplified approach 
Poisson's ratio 


stress function 
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Subscripts: 


A quantities associated with prebuckling state 

B quantities associated with buckling state 

BUCKLING OF THE EQUATORIAL SEGMENT 


Buckling Equations and Boundary Conditions 

A narrow equatorial segment of a thin-walled spherical shell is shown in figure 2. 
This shell segment is simply supported by vertical line reactions at its upper and lower 
edges. It is of interest to determine how thick this spherical segment must be to support 
its own weight without buckling. 

For the equatorial segment the nonlinear equilibrium equations which include body- 
force terms can be written as follows: 


9N X 8N*y 

9x 8y 


+ Px = 0 


(la) 



(lb) 


D v% + Nx + N y 

R 


-L(n x M 

8x\ x 8xJ 8y ' 



" f Xy W +Pn = ° 


(lc) 

The coordinate system of the shell segment is shown in figure 3; the y coordinate is the 
equator of the spherical shell. Membrane stress resultants, displacements, and body 
forces are also indicated in figure 3. Equations (1) are equivalent to equations (6) of ref- 
erence 1 except that body forces are not included in reference 1. To obtain buckling equa- 
tions, the stress resultants and normal displacement in equations (1) are regarded as 
being composed of a prebuckling component and a buckling component, a small change 
which occurs during buckling. Thus, the assumption is made that 


N x = N Xj a + N x> b 

(2a) 

N y = N y> A + N Yj b 

(2b) 

r xy = N X y,A + Nxy ,B 

(2c) 



W = W A + W B 


(2d) 


When equations (2) are substituted into equations (1) and products of buckling quantities 
are disregarded, the following equations are obtained: 


8 N x ,A 

1 — 4- 

8 ^xy,A 

8 N x ,B 
— + 

8N xy,B 

9x 

0 y + 

Ox 

8 y 


8N xy,A 

^B + 

8N xy,B 

9 y 

ax 

0 y 

8x 


+ p x = 0 

(3a) 

+ Py = 0 

(3b) 


Dv4 ( w a + w B ) + Nx,A + Ny,A * Nx,B + Ny> - - - JL (n : 


8W / 




te \ Nxy ’ B 


tel Nxy > A 


CD 

a 

8 y J 

8 y 

awA 

a 

0 y/ 

~ 8 y 

aw B \ 

a 

9 y / 

0 y 


R 

aw/ 


x,A ' 


0W/ 


0x J 8y 
' 8W A N 


- — N v A 

Atr \ y ’ A 


aw^ 
9 y. 


8x / 8x 


ax / ay 


N y,B 


9 y 


N- 


xy,B 


8Wa 
8x / 


aw 


B 


8x 


a f 9W B 
a x \ X ’ A a x / 


+ p n = 0 


a L 9 w b 

9y \ y ’ A 8y j 


(3c) 


Buckling equations can be obtained from equations (3) on the basis of the following 
considerations: (1) The prebuckling terms alone must satisfy the nonlinear equilibrium 

equations (eqs. (1)). (2) As is conventional in shell buckling theory, the prebuckling state 

is assumed to consist of stresses given by the linear membrane theory and normal dis- 
placements which are independent of x and y. The resulting buckling equations are 


8N x,B [ 8N xy,B Q 
8x 8y 


(4a) 


0N 


y,b 


dv4w b , N *' b * N 1^p x ^ 
a R x ax 


B 


9 y 

aw 


, 8N *y,s _ 0 
8x 


(4b) 


+ P. 


y 9y 


- N- 


x,A ' 


8x 2 


N. 


y,A- 


9 \ 

0y 2 


- 2N- 


xy,A 


a 2 w B 

8x 8y 


= 0 


(4c) 
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Equations (4a) and (4b) can be satisfied identically by defining a stress function 4> B 
such that 




Nx ’ B_ a 2 

9 y 


B 


N. 


8 2 4 > 


y,B = 


B 


8x 2 


N- 


'xy,B - “ 


8 ^B 
8x 8y 


The stress function $ B is determined by a compatibility equation obtained from 
reference 3 


v 4 $ b = v 2 w b 


(5) 


The boundary conditions for buckling at the upper and lower edges of this equatorial 
segment are taken to be 


N x ,B = o 

(6a) 

v B = o 

(6b) 

W B = 0 

(6c) 

8 2 W b 
a = 0 
8x 2 

(6d) 

t Nx B + Ny b 

at x = ±£. With these boundary conditions assumed, the term 5 in equa- 

2 fv 

Et 

tion (4c) can be replaced by the quantity — W B , as shown in 

R 2 

appendix A. Since the 

equatorial segment is loaded axisymmetrically and no lateral 
shown, as in appendix B, that 

load is applied, it can be 

II 

1 

* 

(7a) 

Py = ° 

(7b) 
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(7c) 



(7d) 


Nxy,A - 0 (7e) 

where y is the weight per unit volume of the shell-wall material and k is the propor- 
tion of the total shell weight supported at the upper edge. Equation (4c) now becomes 


8W t 


dv 4 w b + “1 w B - yt ax 

R 


- yt 


x + 



+ yt 


2x + 



( 8 ) 


The governing buckling equations are equations (5) and (8) with unknowns the stress 
function 4> B and the normal displacement function W B and with boundary conditions 
given by equations (6). The problem is to determine the characteristic values of the shell 
thickness t for which solutions to equations (5) and (8) exist that satisfy the boundary 
conditions, equations (6). 


Solution of Equations by Galerkin Method 
Truncated series expressions for W B and $ B are assumed as follows: 


M 

w B = c° s a 2 A m sin 


m7T(X + -ii 
2 


m=l 


M 


*B = 


-cos 


ny V 

R Z/ 


Et A 

lt Am 


m=l / n_\^ /m7r\ 2 

w U / 


sin - 


mir|x + |i| 


(9a) 


(9b) 


Equations (9) satisfy equation (5) and the boundary conditions given by equations (6) 
identically. They also satisfy equation (8) identically in the y-direction. 
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Now, by the Galerkin method equation (10) is multiplied by 



(q = 1,2,3,. . .,M) 


and the resulting equation is integrated over the interval 


- — = x = — to yield 

Li Li 
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For determining the critical-thickness parameter of a shell, it is convenient to 
express equation (11) in terms of nondimensional parameters as follows: 


" X ( Ac L 


f“( k - i) 

n 2 - (q7T/3) 2 - 1 

l & 2 

+ (cltt/3)?) 2 

M 

2n 2 - ( mn/3 ) 

_ y 

4a 

L 

m=l 

I"n2 


- P 


i2(q2 - m^) qm 


m+q odd 
where 


[n 2 + (q7r/3) 2 ] 


q2 - m 2 


Am - 0 


y L 

Ot = - 7 — 

E 




(q= 1,2,3,. . .,M) (12) 


X = 


12(l - p 2 )r 2 


Equation (12) represents M homogeneous algebraic equations. In matrix notation, 
these equations are 


[ a qm] " X1 {A} = 0 


are 


where the elements of j^qrnj 

«(k - n2 - (q^/3) 2 ] - 1 


q m 


a qm — 4(2- 


a qm - 0 


[n 2 + (q7r/3) 2 ] 
2n 2 - (m7r|3)^| 


2 ] 2 _ 

Mq 2 - 


(q 2 - m 2 ) 


- p 


(13) 


(m = q) 


qm (m q 


[n 2 + (q7r/3) 2 ]' 


q 2 - m 2 \ m + q odd J 


\m + q even/ 
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and 


(Z ^ 
A 1 
A 2 


<» = 


Now, by noting that 

2n 2 - (m7rj3) 2 


a M 


2 j3 2 = 4n 2 - (m 7 TjS) 2 - (qrr/3) 2 


7r 2 (q 2 - m 2 ) 

equation (13) can be put in the more convenient symmetric form 

[bqm] - XI © = 0 

where the elements of jbqmj 


?r 2 (q 2 - m 2 ) 


are 


“ |k -f 


bqm = bmq = 


bqm - bjnq — 4a 


bqm - b m q - 0 


n 2 - (q tt/3) 2 


- 1 


[n 2 + (q7r/3) 2 J 

4n2 - (m77/3 

2 

) 2 - (q^) 2 

qm 

n2 + (m7T/3) 2 

n2 + (q7r/3) 2 

7T 2 1 

- mf 


(m = q) 


fm ^ q ) 

1 m + q odd 


m/q ] 

l m + q even 


(14) 
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and 



Equation (14) is a standard eigenvalue problem, the nontrivial solutions of which 
yield M characteristic values of the thickness parameter X as a function of a, /3, 
k, and n. For a given a, and k, the limiting value of X obtained by maximizing 
with respect to n as M — yields the critical value X cr ; that is, for X greater 
than X cr the shell does not buckle. In practice this critical value of X can be deter- 
mined by increasing M until the maximum value of X with respect to n converges 
to its limiting value. 

NUMERICAL RESULTS AND DISCUSSION 

Critical values of the thickness parameter X have been determined over the range 
of the specific -weight parameter 10"^ = a = 10“^ for radius-length ratios /3 of 3, 5, 
and 10 and for upper-edge support conditions of k = 0, 1/4, 1/2, or 1. Since most metal- 
lic materials have a y/E of the order 10“® per inch (4 x 10“^ per meter), and most 
plastics have a y/E of the order 10" ^ per inch (4 x 10~6 per meter), a value of a 
equal to lO"? might indicate a metallic shell 10 inches (25 cm) long, and a value of a 
equal to 10“4 might indicate a plastic shell 1000 inches (25 meters) long. 

The number of terms in the truncated series solution M required for convergence 
varied from 30 for a = 10~4 to 70 for a = 10“^. For support conditions where k = 1 
or 1/2, convergence was found to be reasonably rapid. For k = 1 and a = 10“^, typical 
values of n at which X became maximum were 169 for /3 = 10 and 137 for /3 = 3; 
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for k = 1/2 and a = 10 -4 , X became maximum at n = 213 for /3 = 10 and at 
n = 170 for /3 = 3. 

For the conditions where k = 1/4 or 0, however, convergence was found to be 
extremely slow when the truncated series solution was taken in the form of equations (9). 
This difficulty was removed by using a modified truncated series solution of the form 


N+M 


W B = cos ^ X Am Sln 


m7r|x +-|i 


m=N 


The eigenvalue X was then maximized with respect to both n and N. For k = 1/4 
or 0, this eigenvalue maximized at n = 0 while the value of N varied from 1 for 
a = 10 -4 to 518 for a = 10“7, and the value of M again varied from 30 to 70. 

When more than one -quarter of the weight of the equatorial segment is supported 
at the upper edge of the segment, the buckling mode tends to be composed of few waves 
in the axial direction (x-direction) and many waves circumferentially (y-direction). When 
more than three-quarters of the segment weight is supported at the lower edge, on the 
other hand, the buckling mode tends to be axisymmetric (no waves circumferentially) and 
composed of a few or many waves in the axial direction. 


The critical-thickness parameter obtained by the simplified approach X* is 
determined in appendix B as 


x *-(h) 2< * 


(OSk^ 1/4) 
(1/4 SkSl) 


The ratio X cr /x* as a function of the specific-weight parameter a for radius -length 
ratios /3 of 3, 5, and 10 is presented in figures 4, 5, 6, and 7 for k = 0, 1/4, 1/2, and 1, 
respectively. The ratios X*/af2 and X cr /a2 are presented as functions of k for 
a = 10-7, 10"®, 10 -5 , and 10" 4 in figures 8, 9, and 10 for j3 = 3, 5, and 10, respectively. 

The ratio X cr /x* is less than unity for all conditions considered in this report as 
shown in figures 4, 5, 6, and 7. These figures also show that ^-cr/^* increases as the 
radius-length ratio and the specific -weight parameter decrease. In addition, these fig- 
ures show that X cr/ /x* is significantly lower for k = 1 or 1/2 than for k = 1/4 or 0. 
The proportion of the shell weight which should be supported at the upper edge to mini- 
mize the shell thickness (or the tendency to buckle) is indicated by the lowest points on 
the curves for X*/a^ and *-cr/°? in figures 8, 9, and 10. The lowest value of \*/a% 
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is always at k = 1/4 regardless of the value of a whereas the lowest values of 
Xcr/^ always occur at values of k greater than 1/4. The optimum value of k for 
Xcr/* 2 increases as the radius-length ratio and specific-weight parameter increase. 

For /3 = 10 and a - 10“^ the optimum value of k is a little less than 1/2. 

SUMMARY OF RESULTS 

Linear buckling equations have been used to study the problem of buckling of a 
simply supported equatorial segment of a thin-walled spherical shell loaded in the axial 
direction by its own weight. By using these equations the following results, which can be 
used in the design of an equatorial or rectangular segment of a spherical shell, were 
obtained: 

1. The simplified approach, in which the shell segment is assumed to buckle if the 
maximum compressive stress is greater than the critical compressive stress for a com- 
plete sphere loaded by uniform external pressure, leads to larger critical thicknesses 
than those obtained by the method of this report. Differences between the results of the 
simplified approach and those of the present method are greater for the higher values of 
radius-length ratio and specific-weight parameter; that is, the simplified approach is 
more conservative for these cases. 

2. The optimum support conditions of the shell segment depend on the geometry 
and material properties of the shell. For the range of parameters considered, the opti- 
mum proportion of the shell weight supported at the upper edge increased from approxi- 
mately one-fourth to one-half as the radius-length ratio increased from 3 to 10 and the 
specific-weight parameter increased from 10"^ to 10“^. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., August 29, 1968, 

124-08-06-13-23. 
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APPENDIX A 


MODIFICATION OF THE EQUATION FOR THE EQUILIBRIUM OF 
FORCES NORMAL TO THE SHELL MIDDLE SURFACE 

Equation (4c) represents the equilibrium of forces in the direction normal to the 
middle surface of the shell segment. In this appendix it is shown that with the boundary 

Nv R "I" N v T3 

conditions assumed in the present report, the term — 2 — — — ^ — in equation (4c) can be 

R 

Et 

replaced by the quantity — Wg. 

R 2 

N x b + Ny b 

The relationship between — - — — - — and the stress function $g can be 


R 


written as follows: 


N x ,B + N y , B _ 1 /9 2 4 >b + ^B\ = v2# B 


R 


R \ 0y2 


dx* 


R 


(Al) 


The stress function 4>g is related to the normal displacement Wg through the com- 
patibility equation, equation (5), 

v 4 *B=f v2w B 


or 


v2 v2$ b - f w bJ =° 


(A2) 


Now, for any function \p in a two-dimensional region S, if = 0 over S and 

9 Et 

ip = 0 on all boundaries of S, then p = 0 in S. If V^$g - — Wg can be shown to be 

zero on the boundaries of the shell, then this quantity is also zero throughout the shell 

and Nx?B + Ny > B can be replaced by ^ Wg in the equilibrium equation. 

R 


f WB ^ x = ± i 

lowing stress-resultant-displacement equations obtained from reference 1: 


The value of V^^g - — Wg at x = can be examined with the aid of the fol- 

R 2 
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APPENDIX A 


n x,B 


n y,b 8 

N xy,B 


Et 

9U r 

° i n 

av B 

CM 

1 

T—l 

9x ' M 

8 y 

Et 

0 Vg 

9Ug 

1 - M 2 

_3y 

9x 

Et 


0 Vg^ 


W R i 

+ (i + 


R 


2(1 + m) \ 9 y 9x 


(A3 a) 


(A3b) 


(A3c) 


With the boundary conditions assumed herein (eqs. (6)), N Xj g = Vg = Wg = 0 at x = 

B Vg = 0 at x = ±-t then ^ = 0 at x = ±-^ and from equation (A3a), = 0 

s oy i 9x 

at X = ± J' Furthermore > from equation (A3b), Ny jB = 0 at x = ±j; therefore, 




9 2 $ 


9x £ 


9 y' 


r - f w B = N y ,g + n X)B - f w B = o 


at 



and the conditions for replacement are fulfilled. 
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APPENDIX B 


STRESS DISTRIBUTION DETERMINED BY UNEAR MEMBRANE THEORY 


The linear membrane equilibrium equations for the symmetrically loaded equatorial 
segment shown in figure 2 are solved to determine the prebuckled stress state. When 
N X y = 0 and derivatives with respect to y are zero, these equations are 


9N- 


x,A 


N. 


9x 


x,A + N y,A 
R 


+ Px = 0 


- P n = 0 


(Bl) 


(B2) 


The body-force functions p x and p n are obtained by projecting the external 
loads, herein the weight of the shell, in directions tangent to and normal to the shell. By 

1/x' 2 


assuming J is negligible, the following expressions for p x and p n are obtained: 

P x = -yt (B3) 

(B4) 

Substituting equation (B3) into equation (Bl) and integrating yields 


yt 

Pn= "R X 


n x ,a = ytx + c 

and substituting equations (B4) and (B5) into equation (B2) yields 

hTy,A = -2ytx - C 


(B5) 


(B6) 


where C is a constant determined by the value of N x a at the upper or lower edge of 


1/L 


the shell. By assuming J is negligible, the values of N Xj a at the upper and 
lower edges become 


n x ,a = ytLk 

n x A = ytL(k - 1) 


L’ 

X = — 
2 


x = -w 


(B7) 

(B8) 



APPENDIX B 


Substituting either of these two values for N X; a into equation (B5) and solving for C 
gives 

C = ytL^k - (B9) 

Substituting equation (B9) into equations (B5) and (B6) yields the required membrane- 
stress resultants 


n x ,a = rt 


X + [k 


-i 1 


Ny,A = -yt 


2x + 



(BIO) 

(Bll) 


The maximum compressive-stress resultants are: 
For x = - jf, 


N X ,A = rtL(k - l) 


For x 


L 

2’ 


N y>A --ytLfk + iJ 


(O^kS 1/4) 


(1/4 ik^l) 


By assuming that the spherical shell segment buckles if the maximum compressive stress 
is greater than the critical compressive stress for a complete sphere loaded by uniform 
external pressure, that is, 

■p4.2 

n x,A or N y,A = ’ 

R\|3(l - M 2 ) 

the following critical values for X are obtained: 

x * = (l - i ) 2 " 2 
x *=(M ) 2<t2 


T— 1 

VII 

VII 

o 

(B12) 

(1/4 £ k S 1) 

(B13) 
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Figure 1.- Equatorial segment of spherical shell showing rectangular segment which simulates lunar surface in the 
lunar orbit and landing approach simulator at the Langley Research Center. 
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Figure 2.- Geometry and support conditions of an equatorial segment of a spherical shell. 
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